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Abstract 

We analyze the impact of two equal billiard balls in three ideal situations: when 
the balls freely slide on the plane of the billiard, when they roll without sliding 
and when one of them freely slides and the other rolls. In all the cases we 
suppose that the contact between the balls is smooth. We base our analysis 
on some recent general theoretical results on ideal impacts obtained by means 
of Differential Geometric Impulsive Mechanics. We use symbolic computation 
software to solve the computational difficulties arising by the high number of 
degrees of freedom of the system. Some particular but significative impacts, 
with opportunely assigned left velocities and positions of the balls, are analyzed 
in details. The results admit easy interpretations that turn out to be in good 
agreement with the reasonable forecasts and the behaviours of real systems. 
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Introduction 

The impact problems arising in the study of billiard games has always represented, 
for their possibility of being modelled with a lower or higher degree of simplification 
and the facility of comparison of the theoretical results with the real behaviour of 
the balls in the game, interesting challenges and severe tests for all who study Classi- 
cal Impulsive Mechanics. The geometric approach to these mechanical problems has 
recently given important theoretical results, when techniques of modern Differential 
Geometry were fruitfully applied to investigate the concept of ideality of unilateral 
constraints ([!]), even in presence of kinetic constraints ([2]). Unfortunately, because 
of computational difficulties due to the number of degrees of freedom of the systems, 
often theoretical results cannot be straightly applied to analyze practical problems: 
the impact between two equal billiard balls can be modelled with an impulsive me- 
chanical system with 10 degrees of freedom (two for the position of the center of mass 
and three for the Euler angles for each ball) simultaneously subject to an instanta- 
neous unilateral constraint (the contact condition), possibly together with permanent 
and/or instantaneous kinetic constraints (the rolling conditions), all of them involving 
several parameters. The calculations arising in the study of this impact are almost 
prohibitive if approached with hand techniques. Luckily, nowadays the use of symbolic 
computation softwares can give a very great help to solve (part of) these problems. 

In this paper we present the application of the theoretical results of [3 [2] to the 
classical example of impact of identical billiard balls in three ideal but nevertheless 



significative cases: the smooth case, when both the bahs shde on the horizontal plane; 
the rolling case, when both the balls roll without sliding on the horizontal plane; a 
mixed case, when a rolling ball impacts with a ball that slides. In everyone of these 
cases, the theory is powerful enough to allow a very general approach, with both the 
balls moving on the plane, and with general orientations of the balls at the impact. 
Even in this general situation, some results (for example concerning the angle formed 
by the directions of the balls after the impact) will appear in a very clear form. 
However, in order to have a clearer interpretation of some results, it will be almost 
unavoidable the analysis of particular impacts. This choice will appear even more 
necessary in the mixed case, if we want to avoid almost meaningless situations. In 
these cases, the results assume very simple aspects, that can be easily compared with 
those of well known experiences. 

In the first and simplest case, we had to solve a 11 x 11 linear system (fortunately 
with a very simple structure), while in the second and more complicated case we 
had to solve a 15 x 15 linear system, both involving several parameters, such as the 
trigonometric functions of the Euler angles of the balls: to this aim, we used CoCoA 
4.7, a freeware software for exact symbolic computation ([5]). 

We choose to focus our attention to the most meaningful results from the point 
of view of the physical behaviour of the balls. Then the computational aspect of the 
analysis is not treated in details in the paper. Moreover, we point out that, once 
the results of [H [2] are assumed, the algorithm to determine the behaviour of the 
balls after the impact starting by the knowledge of positions and motions of the balls 
before the impact (is based on but) do not involve refined geometric techniques. Then, 
for the same reason written above, we describe in details the algorithm (that is, we 
describe in details the results of (TJ [2] ) , and then we analyze in details the results 
obtained, but we expose only a very concise description of the geometry underlying 
the problem. 

1 Preliminaries 

The theoretical foundation of ideality for unilateral constraints acting on a mechanical 
system with a finite number of degrees of freedom can be fruitfully framed in the 
context of Modern Differential Geometry, in particular using the concepts of fibre 
bundles, their jet-extensions, subbundles and tangent spaces, additionally endowed 
with suitable metrics. For the reader interested in the general theoretical approach, 
we refer to [1] [2] and the references therein. 

1.1 The geometric environment of the problem 

To study the particular case of two equal billiard balls of radius r and mass m moving 
on the billiard plane, we introduce the configuration space-time V, a differentiable 
manifold fibered over the time line and parameterized by 11 coordinates. We can 
choose these coordinates as (t , , yi , i/'i , "(^i , V^i 1 , 2/2 , ■'/'2 , ^^2 , </72 ) , where t is the time 
coordinate, x^, yi are the coordinates of the center of mass of the i-th ball and ipij'di, ipi 
are the Euler angles giving the orientation of the i-th ball. The simultaneous presence 
of the balls on the plane implies that not all the configurations of V are admissible, 
since at any instant the coordinates must satisfy the relation {x2 — xi)^ + (?/2 ~ 2/i)^ > 
4r^. There is then a natural positional constraint S C V expressed by the condition 

S := {{X2 - xi)2 + (2/2 - yi)' - = 0} (1) 

that must be satisfied in the moment of the impact (and then called the impact 
constraint). 

The absolute-velocity space-time Ji(V) is the first jet-extension of V, parameter- 
ized by 21 coordinates: the first 11 are the coordinates running on V, the remaining 
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10 can be grouped in a vector p = (1, i i, y i, -0 i, ?? i, 99 1, i 2, 2/ 2, '0 2, 2, 95 2): the 
presence of the first constant coordinate stresses the necessity of respecting the affine 
nature of the velocity space-time. 

The space V{V) describing the possible impulses acting on the system is the verti- 
cal fiber bundle of V. ft is a differentiable manifold parameterized by 21 coordinates, 
whose first 11 are the coordinates running on V and the remaining 10 can grouped 
in a vector I = {0, I^^, ly^, I^^, I^^, I^^, I^^ j ^2 j -^lAa j ^1^2)'- t^^is t™e the presence 
of the first constant coordinate is due to the necessity of respecting the vector nature 
of the impulse space, and points out that the fibers of ViV) are the vector spaces 
modelling the affine spaces given by the fibers of Ji(V). 

The relation between the space of velocities and the space of impulses is cleared 
by the very nature of the impact problem: given an input velocity (the so called 
"left-velocity" ) of the system and an impulse I, the sum 



Pi? PL 



I 



is required to be an output velocity ( "right- velocity" ) of the system. A similar 
situation is perfectly framed in the relations between affine spaces and modelling 
vector spaces. However, once the situation is cleared up of possible confusions, 
from the computational point of view the significative parts of velocities p and im- 
pulses I are given by the 10 components (x i,y ijip i,ip i,x 2,il 2,'>P 2,'^ 2,^ 2) and 
, lyi, I^i, I^di, Itpi, Ix2T Iy2T '^i'2! -^•&2T -^^2) respectively. These 10 uples will be used 
without ulterior specifications in the calculations. 

A scalar product $ acting on the space of impulses V{V) is defined by using the 
positive definite matrix g that gives the kinetic energy of the system as a quadratic 
form in the velocity coordinates {x i,y i,^' i,^ i^x^^y 2,ip 2,"^ 2,^ 2)- The ele- 
ments of g are in general functions of the positional coordinates (xi, V'l, (pi, 2:2, y2,^2,'&2,V2) 
and take automatically into account the massive properties of the system. Consid- 
ering the significative part of the impulse vectors and denoting with A the inertia 
momentum of the balls with respect to the center of mass, we have 

\ 



(2) 



/ 

A simple calculation shows that this scalar product is positive definite if and only if 
sindi ^ and sin?92 7^ 0. 

Some subspaces of the velocity space-time JiiV) can be naturally introduced. 
The presence of the constraint ^ determines the space T of the velocities that are 
tangent to the positional constraint. Straightforward arguments (of Differential or 
even Elementary Geometry, since the relative velocities of the centers of mass must 
be orthogonal to the radii of the balls at the instant of the impact) show that in our 
case T is given by those velocities obeying the condition 
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{x 2 ~ X i){x2 - xi) + [y 2 " y i){y2 ~ yi) = 0. 
Moreover, the possible presence of rolling conditions 

X i — r'd i sin -ipi + rip i sin di cos -ipi — 

y i + rd i cos ipi + ripi sin 'di sin ij^i ~ 



(3) 



(4) 
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for one or both the baUs restricts the space of admissible vefocities to a subspace 
A C Ji(V). Note that, when the system is subject to the space A is restricted 
before, during and after the impact. However, due to the hnearity of (jl]), ^ stiU owns 
the geometric structure of affine space. 

2 Ideality of the unilateral constraint 

The geometric setup synthetically described above allows the introduction of a sub- 
space ADT C A Q Ji(V) of the space of admissible velocities A of the system before 
and after the impact given by the admissible velocities of the system when the balls 
are in contact. We assume that the contact between the balls is smooth (otherwise an 
ulterior restriction of the space of admissible velocities of the system when the balls 
are in contact must be introduced for the presence of an ulterior kinetic condition) . 

Due to the affine structure of these spaces and the existence of the metric structure 
given by for every element p in A, there exists a closest element V{p) in ^ n T. 
The element 7'(p) can be thought of as the component of p tangent to the impact 
constraint, while the difference 

i^{p) = P-rip) (5) 

can be thought of as the component of p orthogonal to the impact constraint. 

The most natural behaviour of the impact constraint is then given by a "reflec- 
tion" . The impact with the constraint leaves invariant the tangent component V{pl) 
of the left velocity and reflects the orthogonal component I'^{pl). Then, for every 
admissible left velocity p^, the corresponding right velocity pu is given by 

Pr^Pl-2I^{Pl). (6) 

It can be shown (see [HEj) that the right velocity p^ given by ([6]) belongs to A if 
Pl A. Moreover, it preserves the kinetic energy in the impact (for every observer 
for which the constraint is at rest, that are the only ones for which the conservation 
of kinetic energy has a clear meaning [4J). This property motivates the ideality of the 
chosen reaction. 

3 The algorithm 

The algorithm to determine the right velocity p^^ starting from a left velocity p l can 
be easily described by the following procedure: 

1) we start from an a priori knowledge of the space of admissible velocities A 
(coinciding with the whole Ji(V) when no rolling constraints act on the balls), 
its subspace AdT (determined by the impulse kinetic constraint [3]) and the 
scalar product O 

2) we have as input an assigned left-velocity p^ G A; 

3) we determine the tangent component V{pl) G AdT and the orthogonal com- 
ponent /-"-(pl); 

4) we determine the right -velocity p^ G A using the relation Pi?, Pl — 21-^ (pl)- 

Clearly, from the operative point of view, the core of the algorithm is point 3): 
it consists in the solution of the Least Square Problem determining the projection 
V{pl) of the assigned left velocity p^ ^ A onto the space AdT, together with the 
condition ^ in the smooth case, together with the conditions ^ and ^ for both 
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i = 1,2 in the rolling case, and together with the conditions Q and (H]) for only one 
i e {1,2} for the mixed case. 

Starting from the 10-uple pi ~ (i f , y f , f , f , f, i ^, y ^, ^Z- i? ^, we 
introduce the function of the 10 unknown {x i,y i, (p i, x 2,y 2,'>P 2,''^ 2, ^ 2) 



2 



+ +a({>^^^{) 

+2AcosV'i (^1 -V'f) (^1 -^k) 



+ m(x2- x^Y + ™ (2)2 - 2/2)^ 
2 / . . \ 2 



(7) 



+a(v-2-^^) +A(§2-^k) +A{^2-vky 

+ 2AC0SV'2 {lh2-i^2) {<p2- ^2) ■ 



Using the Lagrange Multipliers Method, in the smooth case we look (neglecting the 
parameters introduced by the LMM) for a 10-uplc minimizing the function 

Csmooth = Hp - PlIP + A [(±2 - i l)(x2 - Xi) + (7)2 - y l)(y2 - yi)] ■ (8) 

In the rolling case we look (once again neglecting the parameters introduced by the 
LMM) for a 10~uple minimizing the function 

J^roU = ||P-PLP + >^ [{i 2 - X i){x2 - Xi) + {y 2 - y l)iy2 ~ yi)] 

+ /Zii 1 — ri? 1 sin'01 + '"'^ 1 sini^i cos V'l^ 

+ Mi2 ( y 1 + rij 1 cos?Ai + ri^ 1 sini?i sin V'l j (9) 

+ M21 ( i 2 — 2 sin -02 + 2 sin 1^2 cos ^^2 j 

+ M22 (y2 + »'i5'2 cosV'2 + '^'^2 sin?92 sin'02 j • 

In the mixed case (supposing the ball nr. 2 rolling without sliding) we look for a 
10-uple minimizing the function 

Cmixed = IIp-plIP + ^ [{x2 ~ X i){x2 - xi) + {y2 ~ y i){y2 ~ yi)] 

+ A*2i ('±2 - ?'ti'2 sinV'2 + r(p2 sini?2 cos-02) (10) 
+ lJ-22 [y 2 + r-i) 2 cos i}2 + r(p2 sin ^2 sin ^2 



In all cases the minimizing 10-uple gives the required projection 7^(pl), and then 
the main object for the analysis of the impact. 



4 The possible ideal impacts 

In this section we present the results of the computation in the three cases described 
above. To perform this computation, we used CoCoA 4.7, a freeware software for 
polynomial symbolic computation ([3])- For each situation, we present some general 
(in the sense of "depending by the whole set of parameters" ) but significative results, 
and some results obtained for particular impacts (in the sense of "having chosen some 
parameters in a simple but significative way" ) . 

4.1 Case I. The smooth situation 

It can be easily proved that, when the contact between the balls and the plane of 
the billiard is smooth, the impact does not affect the behaviour of the "angle" vari- 
ables (0i,i^i,(/9i,?/'2,i^2jV'2) and the unique variations regard the "linear" variables 
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(i 1, y 1, i 2, 2/ 2) of the velocities of the centers of mass of the balls. In particular, we 
have the following results: 



' f (2/2 - Vi? + i 2 (-^2 -xif + {y^- y{){x2 - Xi){y2 - yi) 
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(11) 
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Note that, while the exit velocity depends on the relative position of the centers of 
mass of the balls at the impact, it does not depend on the orientations of the balls at 
the impact. 

Let us denote with Vf ,Vf , i — 1,2 the "usual" left and right velocities of the 
centers of mass of the balls (so that, for example, the component of V^" along the 
X direction is if ) and with u • w the "usual" scalar product. Then the results (fTTj) 
satisfy the relation 

Yf-V^=Vf.V^ (12) 

In particular (|12p implies that, when the impact happens with one of the two balls at 
rest, independently of the impact angle, the exit velocities of the two centers of mass 
are orthogonal. 

The computation shows also that the reactive impulse —21-^ {pl) has (null com- 
ponents relative to the "angle" variables and) components relative to the "linear" 
variables 



r± r-L _ (^2 -ii)ix2-xi)^ + {y^ -y^)ix2-xi){y2-yi) _ _ , 

^xi — ^X2 — 27-2 ^ ' 

(13) 



T± _ r± _ ix2 ~xf){x2-xi)iy2-yi) + iy^ -yf)iy2-yi? , _ . 
^yi — ^y2 — 



Therefore, each reactive impulse relative to a single ball turns out to be parallel to 
the direction determined by the centers of mass of the balls. This property, together 
with the smoothness of the contact between balls and plane), motivate also why the 
impact does not affect the "spin" of the balls (see {5J. 

Moreover, simple tests show the reasonableness of (ITT|) . When = 0, and the 
impact is direct central, then = 0, Vf = . More generally, when Vf" = 0, 
then Vf turns out to be parallel to the direction determined by the centers of mass 
at the impact. 



4.2 Case II. The pure rolling situation 

It can be easily foreseen that, differently from the smooth case, when both the balls 
roll without sliding on the billiard plane, in general all the 10 component representing 
the velocities of the balls change with the impact. The computation support this 
forecast. In particular, the "linear" velocity variables have the same behaviour in the 
smooth and the rolling cases, so that ([TT|) still hold, together with (fT2l [T3 |l . 

Unfortunately, although explicitly computable, the changes of the "angle" velocity 
variables have complicated expressions, without easy interpretation. For example, we 
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have, 



+ [r^i - 0/2 -yi) cosiJi{y2-yif 



+ 



0/2 - yi)sinV'i - (a; 2 -i;i')cosV'i ix2 - xi){y2 



It is then convenient to restrict our analysis to particular situations, such as when the 
left velocity of one ball is fixed, and possibly null, and fixing the angle between Vj 
and the direction determined by the centers of mass. 

For example, for an impact with the ball nr. 1 at rest, parallel to the x- 
direction and angle between V2 and the direction determined by the centers of mass 
of ^, we have 
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Another simple example is given by the direct central impact with the ball nr. 1 
at rest and V2 parallel to the x-direction. We have 



if 
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COS')/'! COSl?! ±2 

sin??! r ' 
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sm^i ; 



L 



COS "01 X 

sini9i r 
if = 0; 

y§ = 

; R ; r COS 1/12 COS t?2 i ? 

= ^2 ; 

Sm 172 

- 0; 
(^f = 0. 

4.3 Case III. The mixed situation 

The mixed situation, when one of the baUs rolls without sliding and the other slides 
on the billiard plane seems a very artful situation. However, taking into account the 
results presented in [6] about the behaviour of a billiard ball moving with friction on 
the billiard plane, the case of a rolling ball impacting with another ball at rest so 
that the second ball starts to move sliding on the billiard plane is the ideal case that 
better models a "real world" situation. Once again, some simple particular impacts 
give clearer outlooks than general results. For example, let us consider once again an 
impact with the ball nr. 1 at rest (on a smooth part of the billiard plane), the ball nr. 
2 rolling on the billiard plane with V2 parallel to the x-direction and angle between 
V2 and the direction determined by the centers of mass of ^ . If we set A = amr^ , 
we obtain 

2 I 



a + 2 



i^f = 0; 
= 0; 



• R _ a+ 1 . L 
a + 2 



il2 — To •^2 
a + 2 



; D : r 1 cos •d2 f n T i 9 

i^2 = ^2 - — ^ h9^ + ^cOS02 

a + 2 sm V2 \ r 



^2 = ({a + l)d^+ ^ COSV2 ) ; 

a + 2 V r ' 
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^2 = fl^2 - (« + 1) — . 

a + 2 sin'i?2 V r ^ J 

However, in this case the most interesting result is given by the scalar product Vf • 
of the right "linear" velocities: a straightforward calculation shows that 

-2 " (a + 2)2 ^"^^J ^U. 

Therefore, in this situation, the exit velocities of the balls are not orthogonal. Since 
this is, for every billiard player, a very well known "real" fact, the result gives an 
ulterior proof of the good agreement of the result of the mixed model to the real 
situation. 

5 Conclusions 

We used a symbolic computation software to obtain detailed results in the study of 
the behaviour of two equal billiard balls colliding in three ideal situations. The use 
of computation software was almost essential in order to obtain effective results by 
the known theoretical procedure, since the high number of degrees of freedom of the 
system involves calculations that could hardly be approached by hand techniques. 

Although the computation were performed in very general conditions, often the 
results obtained can be interpreted only for particular impacts, such as when one of 
the ball is at rest. Therefore, only a few of the numerous interesting impacts where 
analyzed in details. 

The results in the smooth and the pure rolling situations respect the forecasts 
(easier in the smooth case, more complicated in the other), especially about the angle 
formed by the directions of the balls after the impact, and about the spins of the balls 
after the impact. The analysis of the mixed situation, suggested by similar analyses of 
analogous problems that can be found in literature, although the apparent artfulness 
of the case, give results that are in good agreement to well known real behaviours of 
the billiard balls. 

References 

[1] S. Pasquero. Ideality criterion for unilateral constraints in time-dependent impul- 
sive mechanics. Journ. Math. Phys., 46,n. 11:112904 20pp, 2005. 

[2] S. Pasquero. On the simultaneous presence of unilateral and kinetic constraints in 
time-dependent impulsive mechanics. Journ. Math. Phys., 47,n.8:082903-19pp, 
2006. 

[3] CoCoA Team, CoCoA: a system for doing computations in Commutative Algebra, 
Available at http : //cocoa . dima . unige . it 

[4] S. Pasquero. On Carnot's theorem in time-dependent impulsive mechanics. Extr. 

Math., 1:87-97, 2005. 

[5] A. Salazar and A. Sanchez-Lavega. Motion of a ball on a rough horizontal surface 
after being struck by a tapering rod. Eur. J. Phys., 11:228 232, 1990. 

[6] M. de la Torre Juarez. The effect of impulsive forces on a system with friction: 
the example of billiard game. Eur. J. Phys., 15:184-190, 1994. 



9 



